Diffusion tensor brain imaging at 0.55T: A feasibility study.

PURPOSE: To investigate the feasibility of diffusion tensor brain imaging at 0.55T with comparisons against 3T. METHODS: Diffusion tensor imaging data with 2 mm isotropic resolution was acquired on a cohort of five healthy subjects using both 0.55T and 3T scanners. The signal-to-noise ratio (SNR) of the 0.55T data was improved using a previous SNR-enhancing joint reconstruction method that jointly reconstructs the entire set of diffusion weighted images from k-space using shared-edge constraints. Quantitative diffusion tensor parameters were estimated and compared across field strengths. We also performed a test-retest assessment of repeatability at each field strength. RESULTS: After applying SNR-enhancing joint reconstruction, the diffusion tensor parameters obtained from 0.55T data were strongly correlated ( R 2 ≥ 0 . 70 $$ {R}^2\ge 0.70 $$ ) with those obtained from 3T data. Test-retest analysis showed that SNR-enhancing reconstruction improved the repeatability of the 0.55T diffusion tensor parameters. CONCLUSION: High-resolution in vivo diffusion MRI of the human brain is feasible at 0.55T when appropriate noise-mitigation strategies are applied.

The "hidden noise" problem in MR image reconstruction.

PURPOSE: The performance of modern image reconstruction methods is commonly judged using quantitative error metrics like root mean squared-error and the structural similarity index, which are calculated by comparing reconstructed images against fully sampled reference data. In practice, the reference data will contain noise and is not a true gold standard. In this work, we demonstrate that the "hidden noise" present in reference data can substantially confound standard approaches for ranking different image reconstruction results. METHODS: Using both experimental and simulated k-space data and several different image reconstruction techniques, we examined whether there was correlation between performance metrics obtained with typical noisy reference data versus those obtained with higher-quality reference data. RESULTS: For conventional performance metrics, the reconstructions that matched best with the higher-quality reference data were substantially different from the reconstructions that matched best with typical noisy reference data. This leads to suboptimal reconstruction results if the performance with respect to noisy reference data is used to select which reconstruction methods/parameters to employ. These issues were reduced when employing alternative error metrics that better account for noise. CONCLUSION: Reference data containing hidden noise can substantially mislead the ranking of image reconstruction methods when using conventional error metrics, but this issue can be mitigated with alternative error metrics.

New Theory and Faster Computations for Subspace-Based Sensitivity Map Estimation in Multichannel MRI.

Sensitivity map estimation is important in many multichannel MRI applications. Subspace-based sensitivity map estimation methods like ESPIRiT are popular and perform well, though can be computationally expensive and their theoretical principles can be nontrivial to understand. In the first part of this work, we present a novel theoretical derivation of subspace-based sensitivity map estimation based on a linear-predictability/structured low-rank modeling perspective. This results in an estimation approach that is equivalent to ESPIRiT, but with distinct theory that may be more intuitive for some readers. In the second part of this work, we propose and evaluate a set of computational acceleration approaches (collectively known as PISCO) that can enable substantial improvements in computation time (up to  âˆ¼ 100× in the examples we show) and memory for subspace-based sensitivity map estimation.

On Ambiguity in Linear Inverse Problems: Entrywise Bounds on Nearly Data-Consistent Solutions and Entrywise Condition Numbers.

Ill-posed linear inverse problems appear frequently in various signal processing applications. It can be very useful to have theoretical characterizations that quantify the level of ill-posedness for a given inverse problem and the degree of ambiguity that may exist about its solution. Traditional measures of ill-posedness, such as the condition number of a matrix, provide characterizations that are global in nature. While such characterizations can be powerful, they can also fail to provide full insight into situations where certain entries of the solution vector are more or less ambiguous than others. In this work, we derive novel theoretical lower- and upper-bounds that apply to individual entries of the solution vector, and are valid for all potential solution vectors that are nearly data-consistent. These bounds are agnostic to the noise statistics and the specific method used to solve the inverse problem, and are also shown to be tight. In addition, our results also lead us to introduce an entrywise version of the traditional condition number, which provides a substantially more nuanced characterization of scenarios where certain elements of the solution vector are less sensitive to perturbations than others. Our results are illustrated in an application to magnetic resonance imaging reconstruction, and we include discussions of practical computation methods for large-scale inverse problems, connections between our new theory and the traditional Cramér-Rao bound under statistical modeling assumptions, and potential extensions to cases involving constraints beyond just data-consistency.

Single breath-hold CINE imaging with combined simultaneous multislice and region-optimized virtual coils.

PURPOSE: To investigate the feasibility of combining simultaneous multislice (SMS) and region-optimized virtual coils (ROVir) for single breath-hold CINE imaging. METHOD: ROVir is a recent virtual coil approach that allows reduced-field of view (FOV) imaging by localizing the signal from a region-of-interest (ROI) and/or suppressing the signal from unwanted spatial regions. In this work, ROVir is used for reduced-FOV SMS bSSFP CINE imaging, which enables whole heart CINE with a single breath-hold acquisition. RESULTS: Reduced-FOV CINE with either SMS-only or ROVir-only resulted in significant aliasing, with severely reduced image quality when compared to the full FOV reference CINE, while the visual appearance of aliasing was substantially reduced with the proposed SMS+ROVir. The end diastolic volume, end systolic volume, and ejection fraction obtained using the proposed approach were similar to the clinical reference (correlations of 0.92, 0.94, and 0.88, respectively with p < 0 . 05 $$ p<0.05 $$ in each case, and biases of 0.1, 1.6 mL, and - 0 . 6 % $$ -0.6\% $$ , respectively). No statistically significant differences for these parameters were found with a Wilcoxon rank test (p = 0.96, 0.20, and 0.40, respectively). CONCLUSION: We demonstrated that reduced-FOV CINE imaging with SMS+ROVir enables single breath-hold whole-heart imaging without compromising visual image quality or quantitative cardiac function parameters.

BrainSuite BIDS App: Containerized Workflows for MRI Analysis.

There has been a concerted effort by the neuroimaging community to establish standards for computational methods for data analysis that promote reproducibility and portability. In particular, the Brain Imaging Data Structure (BIDS) specifies a standard for storing imaging data, and the related BIDS App methodology provides a standard for implementing containerized processing environments that include all necessary dependencies to process BIDS datasets using image processing workflows. We present the BrainSuite BIDS App, which encapsulates the core MRI processing functionality of BrainSuite within the BIDS App framework. Specifically, the BrainSuite BIDS App implements a participant-level workflow comprising three pipelines and a corresponding set of group-level analysis workflows for processing the participant-level outputs. The BrainSuite Anatomical Pipeline (BAP) extracts cortical surface models from a T1-weighted (T1w) MRI. It then performs surface-constrained volumetric registration to align the T1w MRI to a labeled anatomical atlas, which is used to delineate anatomical regions of interest in the MRI brain volume and on the cortical surface models. The BrainSuite Diffusion Pipeline (BDP) processes diffusion-weighted imaging (DWI) data, with steps that include coregistering the DWI data to the T1w scan, correcting for geometric image distortion, and fitting diffusion models to the DWI data. The BrainSuite Functional Pipeline (BFP) performs fMRI processing using a combination of FSL, AFNI, and BrainSuite tools. BFP coregisters the fMRI data to the T1w image, then transforms the data to the anatomical atlas space and to the Human Connectome Project's grayordinate space. Each of these outputs can then be processed during group-level analysis. The outputs of BAP and BDP are analyzed using the BrainSuite Statistics in R (bssr) toolbox, which provides functionality for hypothesis testing and statistical modeling. The outputs of BFP can be analyzed using atlas-based or atlas-free statistical methods during group-level processing. These analyses include the application of BrainSync, which synchronizes the time-series data temporally and enables comparison of resting-state or task-based fMRI data across scans. We also present the BrainSuite Dashboard quality control system, which provides a browser-based interface for reviewing the outputs of individual modules of the participant-level pipelines across a study in real-time as they are generated. BrainSuite Dashboard facilitates rapid review of intermediate results, enabling users to identify processing errors and make adjustments to processing parameters if necessary. The comprehensive functionality included in the BrainSuite BIDS App provides a mechanism for rapidly deploying the BrainSuite workflows into new environments to perform large-scale studies. We demonstrate the capabilities of the BrainSuite BIDS App using structural, diffusion, and functional MRI data from the Amsterdam Open MRI Collection's Population Imaging of Psychology dataset.

High-fidelity, high-spatial-resolution diffusion magnetic resonance imaging of ex vivo whole human brain at ultra-high gradient strength with structured low-rank echo-planar imaging ghost correction.

Diffusion magnetic resonance imaging (dMRI) of whole ex vivo human brain specimens enables three-dimensional (3D) mapping of structural connectivity at the mesoscopic scale, providing detailed evaluation of fiber architecture and tissue microstructure at a spatial resolution that is difficult to access in vivo. To account for the short T2 and low diffusivity of fixed tissue, ex vivo dMRI is often acquired using strong diffusion-sensitizing gradients and multishot/segmented 3D echo-planar imaging (EPI) sequences to achieve high spatial resolution. However, the combination of strong diffusion-sensitizing gradients and multishot/segmented EPI readout can result in pronounced ghosting artifacts incurred by nonlinear spatiotemporal variations in the magnetic field produced by eddy currents. Such ghosting artifacts cannot be corrected with conventional correction solutions and pose a significant roadblock to leveraging human MRI scanners with ultrahigh gradients for ex vivo whole-brain dMRI. Here, we show that ghosting-correction approaches that correct for either polarity-related ghosting or shot-to-shot variations in a separate manner are suboptimal for 3D multishot diffusion-weighted EPI experiments in fixed human brain specimens using strong diffusion-sensitizing gradients on the 3-T Connectom MRI scanner, resulting in orientationally biased dMRI estimates. We apply a recently developed advanced k-space reconstruction method based on structured low-rank matrix (SLM) modeling that handles both polarity-related ghosting and shot-to-shot variation simultaneously, to mitigate artifacts in high-angular resolution multishot dMRI data acquired in several fixed human brain specimens at 0.7-0.8-mm isotropic spatial resolution using b-values up to 10,000 s/mm2 and gradient strengths up to 280 mT/m. We demonstrate the improved mapping of diffusion tensor imaging and fiber orientation distribution functions in key neuroanatomical areas distributed across the whole brain using SLM-based EPI ghost correction compared with alternative techniques.

Efficient iterative solutions to complex-valued nonlinear least-squares problems with mixed linear and antilinear operators.

We consider a setting in which it is desired to find an optimal complex vector x ∈ C N that satisfies A (x) ≈ b in a least-squares sense, where b ∈ C M is a data vector (possibly noise-corrupted), and A (·) : C N → C M is a measurement operator. If A (·) were linear, this reduces to the classical linear least-squares problem, which has a well-known analytic solution as well as powerful iterative solution algorithms. However, instead of linear least-squares, this work considers the more complicated scenario where A (·) is nonlinear, but can be represented as the summation and/or composition of some operators that are linear and some operators that are antilinear. Some common nonlinear operations that have this structure include complex conjugation or taking the real-part or imaginary-part of a complex vector. Previous literature has shown that this kind of mixed linear/antilinear least-squares problem can be mapped into a linear least-squares problem by considering x as a vector in R 2N instead of C N . While this approach is valid, the replacement of the original complex-valued optimization problem with a real-valued optimization problem can be complicated to implement, and can also be associated with increased computational complexity. In this work, we describe theory and computational methods that enable mixed linear/antilinear least-squares problems to be solved iteratively using standard linear least-squares tools, while retaining all of the complex-valued structure of the original inverse problem. An illustration is providedtodemonstratethatthisapproachcansimplifytheimplementationandreduce the computational complexity of iterative solution algorithms.

A hybrid high-resolution anatomical MRI atlas with sub-parcellation of cortical gyri using resting fMRI.

We present a new high-quality, single-subject atlas with sub-millimeter voxel resolution, high SNR, and excellent gray-white tissue contrast to resolve fine anatomical details. The atlas is labeled into two parcellation schemes: 1) the anatomical BCI-DNI atlas, which is manually labeled based on known morphological and anatomical features, and 2) the hybrid USCBrain atlas, which incorporates functional information to guide the sub-parcellation of cerebral cortex. In both cases, we provide consistent volumetric and cortical surface-based parcellation and labeling. The intended use of the atlas is as a reference template for structural coregistration and labeling of individual brains. A single-subject T1-weighted image was acquired five times at a resolution of 0.547 mm × 0.547 mm × 0.800 mm and averaged. Images were processed by an expert neuroanatomist using semi-automated methods in BrainSuite to extract the brain, classify tissue-types, and render anatomical surfaces. Sixty-six cortical and 29 noncortical regions were manually labeled to generate the BCI-DNI atlas. The cortical regions were further sub-parcellated into 130 cortical regions based on multi-subject connectivity analysis using resting fMRI (rfMRI) data from the Human Connectome Project (HCP) database to produce the USCBrain atlas. In addition, we provide a delineation between sulcal valleys and gyral crowns, which offer an additional set of 26 sulcal subregions per hemisphere. Lastly, a probabilistic map is provided to give users a quantitative measure of reliability for each gyral subdivision. Utility of the atlas was assessed by computing Adjusted Rand Indices (ARIs) between individual sub-parcellations obtained through structural-only coregistration to the USCBrain atlas and sub-parcellations obtained directly from each subject's resting fMRI data. Both atlas parcellations can be used with the BrainSuite, FreeSurfer, and FSL software packages.

On the shape of convolution kernels in MRI reconstruction: Rectangles versus ellipsoids.

PURPOSE: Many MRI reconstruction methods (including GRAPPA, SPIRiT, ESPIRiT, LORAKS, and convolutional neural network [CNN] methods) involve shift-invariant convolution models. Rectangular convolution kernel shapes are often chosen by default, although ellipsoidal kernel shapes have potentially appealing theoretical characteristics. In this work, we systematically investigate the differences between different kernel shape choices in several contexts. THEORY: It is well-understood that a rectangular region of k-space is associated with anisotropic spatial resolution, while ellipsoidal regions can be associated with more isotropic resolution. Further, for a fixed spatial resolution, ellipsoidal kernels are associated with substantially fewer parameters than rectangular kernels. These characteristics suggest that ellipsoidal kernels may have certain advantages over rectangular kernels. METHODS: We used real retrospectively undersampled k-space data to empirically study the characteristics of rectangular and ellipsoidal kernels in the context of seven methods (GRAPPA, SPIRiT, ESPIRiT, SAKE, LORAKS, AC-LORAKS, and CNN-based reconstructions). RESULTS: Empirical results suggest that both kernel shapes can produce reconstructed images with similar error metrics, although the ellipsoidal shape can often achieve this with reduced computation time and memory usage and/or fewer model parameters. CONCLUSION: Ellipsoidal kernel shapes may offer advantages over rectangular kernel shapes in various MRI applications.

PALMNUT: An Enhanced Proximal Alternating Linearized Minimization Algorithm with Application to Separate Regularization of Magnitude and Phase.

We introduce a new algorithm for complex image reconstruction with separate regularization of the image magnitude and phase. This optimization problem is interesting in many different image reconstruction contexts, although is nonconvex and can be difficult to solve. In this work, we first describe a novel implementation of the previous proximal alternating linearized minimization (PALM) algorithm to solve this optimization problem. We then make enhancements to PALM, leading to a new algorithm named PALMNUT that combines the PALM together with Nesterov's momentum and a novel approach that relies on uncoupled coordinatewise step sizes derived from coordinatewise Lipschitz-like bounds. Theoretically, we establish that a version of PALMNUT (without Nesterov's momentum) monotonically decreases the objective function, guaranteeing convergence of the cost function value. Empirical results obtained in the context of magnetic resonance imaging demonstrate that PALMNUT has computational advantages over common existing approaches like alternating minimization. Although our focus is on the application to separate magnitude and phase regularization, we expect that the same approach may also be useful in other nonconvex optimization problems with similar objective function structure.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects.

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T1 , T2 , and T2∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

Corticospinal Tract Microstructure Predicts Distal Arm Motor Improvements in Chronic Stroke.

BACKGROUND AND PURPOSE: The corticospinal tract (CST) is a crucial brain pathway for distal arm and hand motor control. We aimed to determine whether a diffusion tensor imaging (DTI)-derived CST metric predicts distal upper extremity (UE) motor improvements in chronic stroke survivors. METHODS: We analyzed clinical and neuroimaging data from a randomized controlled rehabilitation trial. Participants completed clinical assessments and neuroimaging at baseline and clinical assessments 4 months later, postintervention. Using univariate linear regression analysis, we determined the linear relationship between the DTI-derived CST fractional anisotropy asymmetry (FAasym) and the percentage of baseline change in log-transformed average Wolf Motor Function Test time for distal items (ΔlnWMFT-distal_%). The least absolute shrinkage and selection operator (LASSO) linear regressions with cross-validation and bootstrapping were used to determine the relative weighting of CST FAasym, other brain metrics, clinical outcomes, and demographics on distal motor improvement. Logistic regression analyses were performed to test whether the CST FAasym can predict clinically significant UE motor improvement. RESULTS: lnWMFT-distal significantly improved at the group level. Baseline CST FAasym explained 26% of the variance in ΔlnWMFT-distal_%. A multivariate LASSO model including baseline CST FAasym, age, and UE Fugl-Meyer explained 39% of the variance in ΔlnWMFT-distal_%. Further, CST FAasym explained more variance in ΔlnWMFT-distal_% than the other significant predictors in the LASSO model. DISCUSSION AND CONCLUSIONS: CST microstructure is a significant predictor of improvement in distal UE motor function in the context of an UE rehabilitation trial in chronic stroke survivors with mild-to-moderate motor impairment.Video Abstract available for more insight from the authors (see the Video, Supplemental Digital Content 1, available at: http://links.lww.com/JNPT/A350).

Local perturbation responses and checkerboard tests: Characterization tools for nonlinear MRI methods.

PURPOSE: Modern methods for MR image reconstruction, denoising, and parameter mapping are becoming increasingly nonlinear, black-box, and at risk of "hallucination." These trends mean that traditional tools for judging confidence in an image (visual quality assessment, point-spread functions (PSFs), g-factor maps, etc.) are less helpful than before. This paper describes and evaluates an approach that can help with assessing confidence in images produced by arbitrary nonlinear methods. THEORY AND METHODS: We propose to characterize nonlinear methods by examining the images they produce before and after applying controlled perturbations to the measured data. This results in functions known as local perturbation responses (LPRs) that can provide useful insight into sensitivity, spatial resolution, and aliasing characteristics. LPRs can be viewed as generalizations of classical PSFs, and are are very flexible-they can be applied to arbitary nonlinear methods and arbitrary datasets across a range of different reconstruction, denoising, and parameter mapping applications. Importantly, LPRs do not require a ground truth image. RESULTS: Impulse-based and checkerboard-pattern LPRs are demonstrated in image reconstruction and denoising scenarios. We observe that these LPRs provide insights into spatial resolution, signal leakage, and aliasing that are not available with other methods. We also observe that popular reference-based image quality metrics (eg, mean-squared error and structural similarity) do not always correlate with good LPR characteristics. CONCLUSIONS: LPRs are a useful tool that can be used to characterize and assess confidence in nonlinear MR methods, and provide insights that are distinct from and complementary to existing quality assessments.

Region-optimized virtual (ROVir) coils: Localization and/or suppression of spatial regions using sensor-domain beamforming.

PURPOSE: In many MRI scenarios, magnetization is often excited from spatial regions that are not of immediate interest. Excitation of uninteresting magnetization can complicate the design of efficient imaging methods, leading to either artifacts or acquisitions that are longer than necessary. While there are many hardware- and sequence-based approaches for suppressing unwanted magnetization, this paper approaches this longstanding problem from a different and complementary angle, using beamforming to suppress signals from unwanted regions without modifying the acquisition hardware or pulse sequence. Unlike existing beamforming approaches, we use a spatially invariant sensor-domain approach that can be applied directly to raw data to facilitate image reconstruction. THEORY AND METHODS: We use beamforming to linearly mix a set of original coils into a set of region-optimized virtual (ROVir) coils. ROVir coils optimize a signal-to-interference ratio metric, are easily calculated using simple generalized eigenvalue decomposition methods, and provide coil compression. RESULTS: ROVir coils were compared against existing coil-compression methods, and were demonstrated to have substantially better signal suppression capabilities. In addition, examples were provided in a variety of different application contexts (including brain MRI, vocal tract MRI, and cardiac MRI; accelerated Cartesian and non-Cartesian imaging; and outer volume suppression) that demonstrate the strong potential of this kind of approach. CONCLUSION: The beamforming-based ROVir framework is simple to implement, has promising capabilities to suppress unwanted MRI signal, and is compatible with and complementary to existing signal suppression methods. We believe that this general approach could prove useful across a wide range of different MRI applications.

Autoregression and Structured Low-Rank Modeling of Sinogram Neighborhoods.

Sinograms are commonly used to represent the raw data from tomographic imaging experiments. Although it is already well-known that sinograms posess some amount of redundancy, in this work, we present novel theory suggesting that sinograms will often possess substantial additional redundancies that have not been explicitly exploited by previous methods. Specifically, we derive that sinograms will often satisfy multiple simple data-dependent autoregression relationships. This kind of autoregressive structure enables missing/degraded sinogram samples to be linearly predicted using a simple shift-invariant linear combination of neighboring samples. Our theory also further implies that if sinogram samples are assembled into a structured Hankel/Toeplitz matrix, then the matrix will be expected to have low-rank characteristics. As a result, sinogram restoration problems can be formulated as structured low-rank matrix recovery problems. Illustrations of this approach are provided using several different (real and simulated) X-ray imaging datasets, including comparisons against a state-of-the-art deep learning approach. Results suggest that structured low-rank matrix methods for sinogram recovery can have comparable performance to state-of-the-art approaches. Although our evaluation focuses on competitive comparisons against other approaches, we believe that autoregressive constraints are actually complementary to existing approaches with strong potential synergies.

Robust autocalibrated structured low-rank EPI ghost correction.

PURPOSE: We propose and evaluate a new structured low-rank method for echo-planar imaging (EPI) ghost correction called Robust Autocalibrated LORAKS (RAC-LORAKS). The method can be used to suppress EPI ghosts arising from the differences between different readout gradient polarities and/or the differences between different shots. It does not require conventional EPI navigator signals, and is robust to imperfect autocalibration data. METHODS: Autocalibrated LORAKS is a previous structured low-rank method for EPI ghost correction that uses GRAPPA-type autocalibration data to enable high-quality ghost correction. This method works well when the autocalibration data are pristine, but performance degrades substantially when the autocalibration information is imperfect. RAC-LORAKS generalizes Autocalibrated LORAKS in two ways. First, it does not completely trust the information from autocalibration data, and instead considers the autocalibration and EPI data simultaneously when estimating low-rank matrix structure. Second, it uses complementary information from the autocalibration data to improve EPI reconstruction in a multi-contrast joint reconstruction framework. RAC-LORAKS is evaluated using simulations and in vivo data, including comparisons to state-of-the-art methods. RESULTS: RAC-LORAKS is demonstrated to have good ghost elimination performance compared to state-of-the-art methods in several complicated EPI acquisition scenarios (including gradient-echo brain imaging, diffusion-encoded brain imaging, and cardiac imaging). CONCLUSIONS: RAC-LORAKS provides effective suppression of EPI ghosts and is robust to imperfect autocalibration data.

Fast submillimeter diffusion MRI using gSlider-SMS and SNR-enhancing joint reconstruction.

PURPOSE: We evaluate a new approach for achieving diffusion MRI data with high spatial resolution, large volume coverage, and fast acquisition speed. THEORY AND METHODS: A recent method called gSlider-SMS enables whole-brain submillimeter diffusion MRI with high signal-to-noise ratio (SNR) efficiency. However, despite the efficient acquisition, the resulting images can still suffer from low SNR due to the small size of the imaging voxels. This work proposes to mitigate the SNR problem by combining gSlider-SMS with a regularized SNR-enhancing reconstruction approach. RESULTS: Illustrative results show that, from gSlider-SMS data acquired over a span of only 25 minutes on a 3T scanner, the proposed method is able to produce 71 MRI images (64 diffusion encoding orientations with b = 1500 s/ mm2 , and 7 images without diffusion weighting) of the entire in vivo human brain with nominal 0.66 mm spatial resolution. Using data acquired from 75 minutes of acquisition as a gold standard reference, we demonstrate that the proposed SNR-enhancement procedure leads to substantial improvements in estimated diffusion parameters compared to conventional gSlider reconstruction. Results also demonstrate that the proposed method has advantages relative to denoising methods based on low-rank matrix modeling. A theoretical analysis of the trade-off between spatial resolution and SNR suggests that the proposed approach has high efficiency. CONCLUSIONS: The combination of gSlider-SMS with advanced regularized reconstruction enables high-resolution quantitative diffusion MRI from a relatively fast acquisition.

Multidimensional correlation spectroscopic imaging of exponential decays: From theoretical principles to in vivo human applications.

Multiexponential modeling of relaxation or diffusion MR signal decays is a popular approach for estimating and spatially mapping different microstructural tissue compartments. While this approach can be quite powerful, it is also limited by the fact that one-dimensional multiexponential modeling is an ill-posed inverse problem with substantial ambiguities. In this article, we present an overview of a recent multidimensional correlation spectroscopic imaging approach to this problem. This approach helps to alleviate ill-posedness by making advantageous use of multidimensional contrast encoding (e.g., 2D diffusion-relaxation encoding or 2D relaxation-relaxation encoding) combined with a regularized spatial-spectral estimation procedure. Theoretical calculations, simulations, and experimental results are used to illustrate the benefits of this approach relative to classical methods. In addition, we demonstrate an initial proof-of-principle application of this kind of approach to in vivo human MRI experiments.